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ABSTRACT 


This thesis deaAs with the digital simulation of 
three phase bridge converter. A detailed representation of 
converter is considered and the mathematical model is developed. 
The generalized equations, describing the state of the converter, 
are derived on the basis of topological considerations. These 
equations can take care of all possible modes of converter 
operation, where the DC link current is continuous. The 
generalized equations are presented in a form that can be 
easily implemented on a digital computer. Computer program 
for converter simulation is developed and is used to analyse 
the disturbances in a tv/o terminal HTOC transmission system. 

In particular, prediction of overcurrents and overvoltages 
caused by the commutation failure and continuous misfire 
in the inverter are given as examples of the application of 
this program. 



OHiPTER 1 


IHTRODUOTIOF 

1.1 GROWTH OP H?DC SISTIHS 

Direct current transmission is one of the exciting 

technical developments in the last two decades. It was 

initiated when the 20 IW Gotland-Sw'eden line was commissioned 

in 1954. With this modest beginning, DC transaission has 

grown in the present decade to an installed capacity of 

more than 2000 HW, in various parts of the world and the 

transmission voltage has increased to +555 SV. There are now 

# 

about 25 projects in operation or under construction through- 
out the world. 

The rapid development of semiconductor technology 
has revolutionised the art of power modulation. Thyristor 
valve is an important component of HVDO system and is used 
in the schen^es that are commissioned after 1970 and has 
replaced the mercury arc valves that v/ere in the previous 
existing schemes. The use of better converter transformers, 
harmonic filters and the powerful control schemes has improved 
the performance of the HVDO systems to a greater extent. The 
development of. HVDO circuit breakers and the introduction of 
multiterminal systems would f\u:ther revolutionise the use of 
DO links in future. 
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1.2 MERIIS MD DmURITS OP HTOG IRIHSMISSIOF 

HVDG transmission technology has developed dne to 
its recognised significant economic and technical advantages 
for certain specific power astern transmission and inter- 
connection applications, lo appreciate this new concept, one should 
make a comparative study of AG and DO transmission. A "brief 
comparative study, bringing out the importance of DC trans- 
mission is presented here in the following text. 

For the same amount of power to he transmitted, a 
bipolar (having two conductors) DO scheme is cheaper than the 
single circuit three-phase AC scheme of similar conductor 
dimensions and clearances. Hence greater power per conductor 
can be transmitted in the case of DG schemes. 

In the case of bipolar DC scheme, the loss of one 
conductor results in a 505^ loss of transmitting capacity 
compared to a complete shut down of transmission in the 
case of a single circuit three phase AG line, thus DC trans- 
mission is more reliable. 

Due to the absence of charging effect, DG cables 
have no limitations on length as against the AO cables. 

It is well known that AG system stability is 
depeident on the power per circuit and the length of the 
line. For long lines it thus leads to the introduction of 
stabilising equipment (series capacitors etc.) which obvlousljr 
increase the systeri cost. In the case of DC, no such 
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stabilising equipment is needed as DC line length has no 
relation to the stability of the system. Moreover, it was 
proved that in the iC/DC systems, the stability of the 
assooiated iC system can be greatly improved through the DC 
pov;er modulation. 

The short circuit level in AO systans may sometime 

demand the upgrading of existing AG circuit breakers where 

as in DO links, the short circuit currents can be limited 
with the help of current controllers. 

The DC link can be used in interconnecting systems 
of differ eit frequencies. Absence of skin effect, less 
corona loss and radio interf erence, and usage of earth as a 
return conductor sire also some of the advantages of DC systems. 

However DC transmission has some disadvantages which 
can be summarised as follows. 

The converters in the DC system generate harmonics and 
thus arising the need for filter circuits. Converters are 
expensive and in addition, filters increaS' the overall cost 
of the system, 

Eventhou^ the DC line itself does not require any 
reactive power eompensation, the converters require much 
reactive power and must be supplied locally. 

Due to lack of S7DC circuit breakers, the tapping of 
DC line is not yet possible and hence the use is limited in 
point to point transmission. 
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It is <|uite clear from the above discussion that the 
line costs are higher for AC and terminal costs are higher 
for DC, So if the line is long enough to give sufficient 
savings in line costs to off-set the higher terminal costs, 
DO provides a more economical solution, 

1.5 imim OF SIMUDAiPIOiT SIIIDISS 

Simulation studies are necessar^r in the planning 
and design stages and in the study of operation of DC links. 
HVDG simulators are used for all these purposes and some of 
the problems that can be simulated with the help of a HVDC 
simulator are as follows i 

1, Ihe development of concepts and equipment for control 
and protection of EVDC system. This inoludesr 

a) Control of power, currait and extinction angle in 
two terminal and multitecminal systems, 

b) The effects of A3 and DC line faults on the control 
perf omance* 

c) Determination of overcurrent aad overvoltage stresses 

in various components of the system. /' 

/ 

2, Aialysis of various phenomena in DCj as well AC/DC 
systems, %vhich includes: 

a) HVDC operation with weak AO systems. 

b) Stability and damping introduced by EHTDO systems. 

c) Control methods for influencing the reactive power 
characteristics of AC systems. 
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d) Smoothing of haimonics through filters and reduction 
of noncharacteristic harmonics. 

So far HVI>C simulator is the most versatile tool 
available for simulaticn studies although digital simulation 
is also being considered for specific problems because of 
its flexibility, A model suffers from the difficulty that 
costs increase as it bee ernes sufficiently versatile to 
represent a wide range of system configurations and also 
scaling of components may produce difficulties, !Ehe dis- 
advantage of the digital simulation for HVDG systems is the 
complexity required in a detailed simulation, where each 
component of the systen has to be modelled accurately, 5?he 
attendant problem of increased computer time and storage has 
so far necessitated in use of different levels of detail in 
the modelling of the system, for example, an analysis of 
dynamic stability of the AO/DC system would require only 
a simplified model of the DO link. But for the transient 
analysis and for accurate simulation of the high frequency 
effects a detailed simulation of the coaverter is essoitial. 

Digital Simulation C onverters : Any system can be simulated 
on digital computer, if the dynamic performance of the qystm 
can be represented by a set of equations, may be algebraic, 
differential or boolean. Digital simulation of bridge 
converters began to be dealt with about fifteen years ago. 

All the simulation raethods so far considered can be classi- 
fied into two groups. 
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7 B 

The first category of simulatioJi methods ^ use a 
traasf ormation technique, by which the state equations 
corresponding to each of the state of the converter are 
formed and solved. The transformation matrix is automatically 
computed by the computer, 

5 6 11 

In the second category^* * , the programs contain 

all the state equations of the possible circuits corresponding 

to the conducting sets of thyristors and at every instant the 

program selects the proper set of eouations for solution. 

In ref. 5, the authors have suggested a method called central 

process method, in which the caitral process starts from the 

instant of firing a valve to the instant of next firing. The 

central processes are mathematically represented by a set of 

differential and boolean equations. Depending upon the state 

of the converter, the differential equations, phase to phase 

voltages are selected and solved, later the work is extended^ 

to simulate the HVDO transmission system under abaormal 

conditions. The central process method seems to be rather 

compli ated vftien compared to the method si^ggested by Htsui 
11 

and Shephei^d , In this simulaticsi work, generalized differ- 
ential equations are formulated after deriving the differential 
equations corresponding to all possible normal combinations 
of conducting valves. 
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1.4 SIWiEr OP IHE THESIS 

la this^the prolslen of converter sinmlation is 

r 

examined with a fresh view point of tying to determine the 
general structure of the system equations for all possible 
modes of operation of the converter including normal and 
abnomal ones. In all, there are 47 possible modes with 
a) two b) three e) four valves conduction. It is shown 
that the system equations for all these modes can be derived 
from a basic set of equations derived from topological 
considerations. Thus the problem of converter simulation 
is put in the general framework of the problem of analysis 
of a network with ^controlled switches. The existence of the 
solution to this general problem is also revealed by exami- 
nation of the basic equations. 

The approach adopted in this thesis is to develop 
sets of generalized differential equations, each set corres- 
ponding to a certain category of the possible modes. This 
approach is along the lines given in ref, 11. However, the 
development of these geaeralized set of equations proceeds 
from the basic equations derived from the network topology. 

A computer program for digital simulation of converters 
in a HYDC link is developed incorporating the features mentioned 
above. This is applied for the transient analysis of the 
converters under both normal and abnormal operating conditions, 
where the direct current is continuous. There are 59 possible 
modes for this condition. 
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Ihe chapterwise description of the thesis is 
given helow; 

i* The second chapter deals mth the mathematical 
representation of the three phase bridge converter and the 
generalized eonations are developed for both normal and 
abnormal operating conditions. The generalized equations 
are presented in a form, that can be easily imiDlemented on 
a digital computer, 

ii. The computer flow chart is included in the third 
chapter and it describes the programming techniques to apply 
the mathanatical model developed in Ghapter 2, for efficient 
canputer simulation. The program is utilized in this chapter 
for the transient analysis of a 5-phase bridge converter, 

iii. The fourth chapter presents the application of the 
ccmputer program for the analysis of abnormal operating condi- 
tions in a two terminal H7D0 system consisting of six pulse 
converters at the rectifier and inverter stations. Two types 
of faults are considered - (i), drop in the inverter voltage 
leaiding to commutation failure and (ii) misoperation of the 
inverter valves. The results of the study are presented and 
ccaapare favourably with those givaa in the previous literature. 

This chapter also presents the developnent of a simple 
procedure for the simulation of converter controls at the 
rectifier. The results with and without the controller are 
compared. 



GHAPTER 2 


modelling oe a coitterier 

2.1 INTRODUCTION 

The varying topology of a converter caused hy the 
cosmencenient and cessation of valve oonduction malces djmamic 
simulation extremely difficult. To simulate dynamic operation, 
the system must he presented in a form acceptable for computer 
solution. Hence to represent the complete system in detail 
by a mathematical model, the e<luations corresponding to all 
possible circuits formed, depending on the conducting states 
of various valves, are considered and the generalised equations 
are formulated. These equations are presented in a form that 
can be easily implemented on a digital computer, 

2.2 EORMDLADION OE EQUATIONS 

The three phase bridge converter is shown in Eig,2.1,a 
in which the valves are numbered in their firing sequence. L^ 
is the smoothing reactor and R^ is its resistance. Eig.2,1.b 
is the network corresponding to the AG power source connected 
to nodes 2, 3 and 4. R*s and L*s are the source resistances 
and inductances respectively. The formulation of equations is 
based on the graph shown in Eig.2,1,c. Here the circuit of 
Eig.2,1,b is. represented by two branches (7 and 9), utiliEiing 
the raultiterminal rexa:esentation for networks. 




Bridqc convei ter cind the associated graph 
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la the formulation of equations valves are treated as 
ideal switches having no forward resistance and infinite 
reverse resistance. The following equations can he easily 
derived from Figs. 2,1,h and 2,1.c. 



The relation between the branch currents and the link currents 
can be expressed in the following form: 


0 0 1 -1 0 0 0 
0 10 0-1 0 0 

-1 0 0 0-1 0 

1 1 0 0 0 0 -1 




The above equation can be expressed in a simple form 


A i , 0 


( 2 . 2 ) 


The tree branch voltages and the link voltages can be related 
as in equation (2,3)* 
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Combining equations (2.1) and (2.3), w© get 





mr 

0 

i (I12P * 

«r 


0 

(%P+\) 

0 


0 

( Ijj^ + I1 2 ) ^ 2 ^ 
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She eijuations (2,2) and (2.4) are the basic equations 
from v/hich the state equations for a given mode can be obtained.. 
In equation (2.2) the valve currents corresponding to the non- 
conducting valves are set equal to zero and the equation can 
be revnritten as 

■%l ^ (2.5) 

where is the reduced matrix obtained from A by eliminating 
the columns corresponding to the nonconducting valves, ij^ is 
the reduced current vector. 

Depending upon the number of conducting valves, 
three cases can be considered; 

(a) Iwo valves conduction: In this case the two valve 
currents and two of the inductor currents can be solved in 
terms of the remaining inductor current, say ig, 

(b) Ihree valves conduction: In this case the three valve 
currents and me of the inductor currents can be solved in 
terms of the remaining two inductor currents, 

(c) Pour valves conduction: In this case all the four valve 
currents can be expressed in terms of all the three inductor 
currents. 

Prom the solution of equation (2.5), one can write 
ij^ w 3? u (2,6) 

where I is the suitable transformation matrix and u is the 
vector of independent inductor currents. 
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Equation (2,4) is reduced tiy putting the voltages corresponding 
to the conducting valves to zero. By deleting the columns 
corresponding to the conducting valves, the matrices and ^2 
are reduced and the equation (2.4) can he reduced to the 
following form; 



In the above equation is the reduced vector formed from 
T after substituting zeros for conducting valve voltages. 

After substituting for ij^ in equation (2.7), we finally get 



It is to be noted that the sum of the dimensions of the vector 
u and Y-g is equal to five and we get five equations in five 
unknowns which have to be solved. Actually what is required 
is, the derivation of the differential equations in u by 
eliminating the valve voltages and algebraic expressions for 
valve voltages in terms of the right hand side quantities of 
equation (2,8), which are known. 

Ihe above procedure is explained in next section by 
considering some normal and abnormal modes of converter 
operation. 
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2.3 mmsms 

Ex, 1 ; A normal mode is considered and let the conducting 
valYes are 1 and 2, let ig he the known current from the 
previous state of the converter. So 

U « ig (2, 9. a) 

As valves 1 and 2 are conducting we can substitute 
i^ = = i^ = ig * 0 in equation (2,2) and = ^2 =* 0 

in equation (2.4). 
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[ (ej+Y^-e^) - (R^+R3+II^) igJ 
Ii, + L^ - 




■* ^\ "* ^i‘:j'^r^’i^^ (Ri +R^+R^ ) I ig 


^6 - ^ \ - V 3 . 


Ex, 2 1 let the conducting valves he 1 , 2 and 3 and let the 
currents ig and are Imown, 

i 

Erom equation (2.5), 
i-j « ig 

^ 2 ’* ^ 8 
H * ^ 8“^9 

Erom equation (2, 5), 1 can he found out and substituting in 

equation (2*6) 

■ % ** 



(2.11) 
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Substituting equaticai (2.11) and the knoMi voltages in equation 
(2.4), wg get 


(b3+lj2)p+(R3+®’2^ ^ 

(li2P+^2 ^ 

0 1 

0 


^8 



(1*^P+R^) 

0 

0 1 

0 


^9 


e 

— (I 2 P+R 2 ) 

(li^ + 112 ) P+(^^-j+R2^ 

0 0 

0 ■ 


‘"'4 ^ 

ss 


0 

0 

1 -1 

0 ! 




0 

0 

0 

0 -1 

1 


^6 




(2.12) 


Prom equation (2.12) the unkncmn vector can be solved and the I 
expressions are as follows: 

pig « D ^ ^ (©^+7^) (lr^+Ii2)“"C®2^1 '*’®1^2^ 1 •* { (R-^+H.^) ^8 

pig = ^ ( 62”‘®1 ^ ^ (R^+B.jj)l*2’“^2^^^^^d^ ^^8 

- {(R^+R2)(L3+%) + R^i.2> igJ . 

^5 * ir '^c 

' R. 

4^ R2l**| *• J i0 (R<j Ii2**R2^'] ) ig 1 « 

^4 ^6 == % 

wli©x*0 »s li^li2 ^2^3 

and D *s + I»^(li^+Ij 2 ^* 
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The above two examples refer to the normal operation of 
converter, in which valve currents are same as the phase 
currents. 2)he following three examples refer to the abnormal 
operation of converter. 

'EXm3 t let the conducting valves be 1 and 4 and ig be the 
known current. 

We can substitute 


as ™ 

and as sa 0 

Prom equations (2.2) and 



( 2 . 6 ) 



in equation (2.2) 
in equation (2.4). 


(2.13) 


Substituting this in equation (2.4), 


^ 0 0 ~1 1 o’ 


'is ■ 



(R^+l^p) 1 _0 1 0 


^2 


e 

0 0-1 0 0 



=2 


o 

> 1 

O 

> 1 

-A 

O 


^5 


0 

0 -1 1 -1 1 I 


.^6_ 




Prom equation (2.14), we get 


■ ■ *** , (o-j^ ' ^'2^'' 


(2.14) 
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“ (e^ - e^) 


-T, 




In this case, it can be noted that there is a short 
circuit on the dc lint and all the three phases are open 
circuited. So the currents through the conducting valves are 
same as the dc link current i.e, , ig« 


Ex.4 t let the conducting valves are 1,4 and 3« Let ig and 
ig be known. 


and 


¥e can substitute 



Prom equation (2,2) 



and also 

+ 1-g 

1-4 * ^8 

ij * ”*^9 


(2.15) 
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item eq.uation (2,7) 


0 


0 

1 

0 


^8 




0 

1 

1 

0 ^ 


ig 


e ■ 

0 

(Xt^ +Ii2)p+(R*j ■^R2) 

_0 

_0 

0 


^2 

BE 


0 

0 

-1 

-1 

0 




! 0 : 

0 

0 

-1 

-1 

1 . 


_'' 6 . 




Prom efiuatioii (2,16), it can "be written as 




( 2 , 16 ) 


^5 * (e2li^+e^Ii2) - e^(L^+l2'^ + ig (R^ Ii2--R2^1 ^ ^ 

* -^5 

Yg « 0 

In this case there is a short circuit on phases 1 and 2 and 
also on dc link, 

Ex«5 t let the conducting valves are 1,2,3 and 4 and let 
ij^ is known. 

So » Y 2 - Tj = 7^ = 0 j ±5 » ig « 0 

and f is a mit matrix. 

Prom equation (2.5) the valve currents can he obtained as 
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. X-| m Xrj 4* i0 ig 
Xn: ** 

- ^ (2.17) 

Prom eguatioa (2.4) 


I 1 0 


Z 


1 0 


0 


i 

t 0 0 

L.- 

t -1 0 
^ -1 1 


-8 


0 

L J 


( 2 . 18 ) 


Prom equation (2,18) the following expressions can be derived: 


1 

pXij sii [ {(e-jl2"^^2^1 ^ (l^*j "^1*2^ ^ ^ ^3 ^^1 ^ ^ 

4- (R*j I i2*^^2^1 ^ 


pis = ir<^o - '‘die) 





[( (6^Ij2+S2^3 ^ 


•“ 6-^ (112+1*^) ) + (5^1j2“^2^3^^7 
. {R.j(l2+I«3)+i^2^3> 


Y5 = Yg » 0 

where D-j* ^-1^2 ^2^3 + ^3^1* 

It car he noted that, in this case, the three phases are short 
circuited and also the do link. 
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2.4 CMSSIFICi£OIOIf OF MODES 

la the normal operation of a three phase "bridge 
converter, either two valves or three valves are conducting 
simultaneously, therefore there will he 12 different modes 
of operation per cycle. In the a"bnormal operation, there are 
27 possible abnormal modes assuming that the DO link is not 
open circuited and all these are considered in the present 
simulation. In the above five examples five different modes 
are explained and in fact all these existing modes whether 
normal or abnormal can be divided into three groups as shown 
in Table 2.1. 


GROUPS 

IT0EM4D MODES 

iSFOEMlD MODES 

TOTilr 

2 Yalves 

( 1 , 2 ) ( 2 , 3 ) ( 3 , 4 ) 

(1,4) (3,6) (5,2) 


conduction 

(4,5) (5,6) (6,1) 

9 

3 Yalves 

(1,2,3) (2,3,4) 

(1,4,3) (1,4,6) 


conduction 

(5,4,5) (4,5,6) 

(1,4,5) (1,4,2) 



(5,6,1) (6,1,2) 

(5.6.1) (5,6,4) 

(5.6.2) (5,6,5) 
(5,2,1) (5,2,4) 
(5,2,5) (5,2,6) 

18 



4 Yalves 


(1,2,5, 4) (2, 5, 4,5) 


conduction 


(5,4, 5,6) (4, 5, 6,1) 

(5. 6.1. 2) (6,1, 2,5) 
(1,5. 5, 4) (1, 5,5,6) 

(1.5. 5. 2) (2,4,6,1) 

(2. 4. 6. 3) (2, 4, 6, 5) 

12 



Table 2,1 1 Glassification of Modes, 
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The generalised equations for the 5 different sets of modes 
are given in the subsequent section in the form that can be 
easily implemented on a digital computer, 

2.5 GENERiilZED BQUillOHS 

2.5»1 G-eneral Definitions! The following points are to be 
noted in connection with the symbols used in the equations, 

i) Bach valve of 5'ig»2, 1,a is represented by ^Im in which 
the first suffix 1, indicates the phase number (1,2, or 5) 
to which the valve is connected and the second suffix m 
indicates whetha: the valve belongs t o upper group or lower 
group. It is to be noted that valves 1,3 and 5 belong to 
upper group and the lower group is formed by valves 2,6 and 4, 

So m can take values 1 or 2 in which 1 corresponds to upper , 
group. 

ii) Valve switching logic: S -,^ » 1 if the valve is conducting 

= 0 if the valve is not 
conducting, 

iii) ill the lower case letters are suffixes except p and e. 
Suffixes except d and c will be defined as per the mode, 

P is a suffix, but p indicates the derivative Cd/dt) 
d is used in connection with smoothing reactor (Ii^) and its 
resistance (R^) and c with dc voltage (V^). 
e is the sinusoidal phase voltage, 

I indicates the current. 
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iv) In the formulation of generalized equations for the 
abnormal cases, the equations for valve currents are to be 
formulated, for which we need the conducting valve numbers. 
For this purpose valve numbers can be obtained from the 
matrix of order 3x2, 



Talve numbers are indicated by n1 , n2, n3 etc., and the above 
matrix elemaats by ¥^2* ^21 » » 

v) 7(l,m) is used to indicate the voltage across a valve, 

where 1 and m are as defined in (i). 


2,5,2 Iv/o Valves Conduction (Formal Mode): 


i. Note i and j for which jDefine p=6-(i+3), 


11. p I a 


v-TT ‘ ’ 


8 3 


iii. Phase currents: * “^S » ** 

iv. V(i,1) * V(j,2) « 0 , 




Is* 




vi. If jsa2 set aa:2j b»1; D=«S=P = -1 and interchange i and 3 
for this step only* 
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If 3362 set a»1j'b*2jD*E = S's1 

(R.+R +R,) 

+ -aq+E^r ^8* 

vii, T(i, 2 ) , 7 ( 3 , 1 ) 

YCP,1>) * V(3,1) - 7(?,a). 


2,5,3 Ihree Valves Oonduction (Formal Mode); 

i. Note t,3 and ^ for "vtoich »: = Sp 2 * 1 and set q=t 

or 


Rote i ,3 and P for which S^ 2 *S ^2 “ ^1 =* ^ I’s-'J. 


In the above step, select i and 3 such that i < 3 , 


ii. If i at 2 get k»i; 1 b3;s = 7;u = 9 
Otherwise set k « 3 ; 1 = i; s = 9; u = 7. 


iii. C >s t ^ 2^3 + ^ 3^1 P = 0 


pig - 5 t{(9ep+Vg)(V^j)-l(ejI.i+ej^I..)} 
- { (Rp+R^)(l^+Il^) + 
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iv. Set B « 1 if p = 2 
= 0 otherwise. 

Set A = -1 if ^ 1 
= 1 otherwise. 

Set X = 3, y = 1, if P = i 
X = 1, y as 3 otherwise. 

’ V ' ^2 = -(Ix * ly’- 

V, If g = 1 set 7(1.1) = 7(3,1) = V(p^ 2) = 0 aad 

S = 2j h sa 1 . 

If a . ~1 set 7(i,2) = 7(3,2) = 7(P,1) = o and 
g » 1; h = 2. 

G 7 ■ I - - 

T(i,g) s= — ^ {aep(L^+ip-(l(e^l^+e^Lp } 

-{ {VRa)(V^3^^®i:h - r^- Is 

- - \h>is 1 


^(Us) » 7(i,g) « 7(P,h). 
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2,5*4 Two Valves Coaduction (Abnormal Mode): 

Let n1 and n2 are the conducting valves, 

I, Note i for which Sii » 

3 « i+1 if i 3 
» 1 if i ^ 5 

P « 6 - (i+3), 

II. pig rs 

as I 2 * = 0 (All the phase currents are zero), 

iii. V(i,1) = V(i,2) = 0 

V(d, 1) « - e^j 

V(3,2) « .V(3,1) 

V(P,1) » - ep 
V(P,2) « «V(P,1). 

iv. Valve curroits: » Ig = 1^2* 

2,5,5 Three Valves Conduction (Abnormal Mode); 

i, JSTote i,3 and ^ for which » S.^ =* Sp 2 * and set q,=1 
or 

Ifote i, 3 and P for which S ^2 “ S ^2 * ® 

In the above, select i aid 3 such that 1 < 3, 
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ii. If i = 2 set s » 7, M = -1, X , 1, y 
otherwise s » 9, M , 1, x = 3, y 

iii. pig » - \ Ig) 

iv. Set A t= 0 if i IS 2 or j B 2 

« 1 othen\?ise 
I_ a» A I 

X -^S 

■ -^6 

I2 » -dx + ly). 


3 

1 - 


V. 


7alve curreats: 

If q ss 1 set a1 

or 


n2 


If q B .t set n1 B jr^2» a2 
Set R B 1 if i „ p 
s= 0 otherwise* 

Ia1 - 8 Ig + q Ij- 


ITji, a5 


ir^2, a3 


JJ. 


'P2 


E 


R1 




-^1 


^a5 “ ^8* 
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Vi, If <i » 1 set V(i,1) , T(o,1) « Y(P,2) = 0 

otherwise set V(i,2) , ¥(3, 2) =. 7(P,1) =0 

If q « 1 and i = P set ¥(3,2) 3 0 otherwise 'V(i,2) 3 0, 

If q»i -1 and i » P set V(3,1) =0 otherwise V(i,1) s 0. 


T(P, 1 ) . -T(P, 2 ), 

# 

2,5.6 Pour Valves Oonductiont 

PI 9 JUj" t { (ejli2+e2^5 ^ ^ (P^l2*"^2^5)^Y 

» {R^ (12+^5) +^2^5 ^ ^ 

■plj * "^C ^ (e^Ii2'’'®2^1 -{R^ ( 1 ^+ 12 ) +1^2^^) 1^ 


+ (R^I.2 - ^2^1> ^9! 

where D ^ l2+Ii2^5+^5^-l » 

ii. Phase ourroits: 1 ^ » -I9 

^2 “ ?9 '*■ ^7 
I3 * 

Set all the voltages across valves to zero. 
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iii, Yalve currents: 

- ’ - 3 3 

Case (a) , j 


l»l 



JTote i, 3 ,k and 1 for which S.j = s « s _ q 

3 I k 2 “ ^2 

subject to conditions i < 3 and k < 1 . 

Sat n1 . , „2 , y ^ 

Ifl.kori,lsetA, 1 ,B = -1,r.3, 

If 3 » k or 3 , 1 set A ^ 0, B „ 1 ^ r , i. 


1 


^n2 “ ^8 • ^n1 

If i » k or 3* = k set 

If i « 1 or 3 » 1 set M * 0 , T = -1 , v * k. 

1^3 - M Ig . I 

^u4 * ^8 ** ^n3* 

Case (b) If , 5 set n1=ir^^, aZ-Sj^, nS-JT,^ and 

note i for which 5^2 ^ set n4«sr^2» 

or 

3 

“ ill ®12 • ^ - *^12' n3 - %2 and 

note if or which « 4 and set n4 » (1*1 

Set a^»1 and other a*s to zero in the following expression s = 

^n1 » ^2 “ 

^n3 " *^^3 ^3% * Ia4 “ I^- 
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2.6 aOFOlUSIOlT 

In this chapter a detailed representation of three 
phase bridge converter is considered. lopological approach 
is usci in fornrulating the generalized equations, ill the 
possible normal and abnormal modes are considered in the 
formulation, 3?he abnormal cases where the dc link current 
is discontinuous are not considered. 



GHiHDER 3 


SBroi>i[DI02T OP COFTERTER 

3,1 IFE-IODUGTIOF 

This chapter presents the digital simulation of a 
3-phase bridge converter based on the mathematical model given 
in Chapter 2. Computer program is developed that utilizes the 
generalized equations for a set of modes. Two sets of equations 
completely describe all the 12 normal modes of operation of a 
converter. 

The program is tested by considering an example of 
transient analysis of a converter feeding into a dc voltage 
source. 

3, 2 3DESGRIPTI0F OE THE COMPUTER PROGRiM 

The computer flow diagram for the simulation is given 
in Pig, 3.1, and is described below in detail. 

i. Input data! Yarious input data that are required are 
mentioned in the flov; diagram, in which initial state of the 
converter corresponds to the numbers of conducting valves at 
time t and initial values of the phase currents, 

ii. Selection of Equations^ Under normal operation of 
converteri two or three valves conduct ^d depending on the 
conducting state of the converter valves the equations that 
are to be solved also change. Prom the formulation of the 




Yes 

/ 


Print the results 
Pig. 3.1 OOKPOTBR PLOW PIAGRM 


Change the state of 
the converter. 
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generalized eluatioas given in previous chapter, the proper 
set of equations is chosen depending upon the converter modes. 

iii. Solution of Selected Equations: Ihe selected differ- 
ential equations are solved using Runge-ICutta method as it is 
the most suitable method for the integration of the equations 
with discontinuities. The accuracy of the solution depends 
on the size of the step and it should he sufficiently small, 
in the simulation a step size of 3 electrical degrees is 
chosen. 

iv, Finding the Discontinuity: Discontinuities are 
created because of tv/o reasons, the commencement and 
cessation of converter valves, 

A valve comes into conducting state only when the 
voltage across it is greater than a certain voltage and a 
firing pulse is present. A valve stops conducting when the 
current through it becomes zero and remains zero for certain 
time. 

Step iii is interrupted if a discontinuity is 
detected. Having detected the presence of discontinuity 
(as per the conditions mentioned above), it is the next 
step to find out its exact occurrence during a given step. 

For this ’linear interpolation* is used as follows. 

Let the discontinuity be because of valve cessation 
i.e. , the current through the valve is becoming negative. 
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^old current through the outgoing valve at time 

previous step and i„^,, he the present 

Uo-U, ii0W 

negative current at time t^ , - 

new- 

Then hy linear interpolation, the exact time t at 
vjhich the current through the outgoing valve becomes zero 
is 

t « t , , + H i -,;,/(i n. - i ) (3.1) 

old old' ^ old new' ^ ' 


whece H is the integration step size* 

Ihe accuracy of this linear interpolation depends 
on the size of the step length chosen. In similar way, 
by using eon, (3* 2), all the dependent variables can he 
interpolated. 


^oia * 


- told> 

(t T, - t ) 

^ old new' 


^^new " ^old^ (5.2) 


V, Updating the State of Converter i After finding out 
the variables at the point of discontinuity , the change in 
the converter state is to be properly implemented for further 
operation. As already maationed, the change may be due to an 
incoming valve or an outgoing valve, implementation of this 
change obviously causes a chaage in the system differential 
equations and these have to be used to perform the next step 
of calculation. 
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3.3 ill EOMPIE 

3.3.1 Description of the SystemJ The example considered 
here is that of a 3-phase bridge converter feeding into a 
voltage source as shown in lig,3,2. It is already assumed 
(in Chapter 2) that the valves have no forv;ard resistance 
and infinite reverse resistance and hence in the simulation 
they are treated as switches 'v^hich are open during the 
blocking state and closed during the conducting state. 

The system parameters are given below; 

Frequency of the ac system 
AC soiirce line to neutral voltage 
Reactance of source inductance 
Source resistance 

Reactance of smoothing reactor(l(^) 

Resistance (%) of smoothing 

reactor 

Value of dc voltage source (Y ) 

C 

Firing angle ( a ) 

All the initial conditions are taken as zero* 

3.3.2 Resultss The results of the siMulation are presented 
in Figs, 5.5 to 3.6 which show 

a) Starting transient of dc link current. 

b) Steady state dc link. current waveform. 

c) Valve voltage waveform in steady state. 


a 50 Hz 

a 1,0 p.U, 
a 0,1 p.U. 
a 0.0 p.U. 

=r 0,1 p.U. 

a 0.005 p.U, 
= 0.65 p.U. 

a 60® 
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d) Rectifier terminal voltage waveform in steady state, 

e) Steady state valve currents in a cycle. 

3.5.3 Discussion: Rig, 3,3 shows the variation of direct fo-veYex^®-. 
current before and after the steady state is reached. It can 
be observed from the graph that the current reached steady 
state in 3 cycles, Ihe average steady state value of the 
current is 5*38 p,u, , whereas the calculated value is 
5.175 P.U., by using equation (3.5) 


dr 


where V 


do 


, 7^^ cosa- C| + Ra)ia 
! ideal no load voltage *= 2,34 


« AO voltage line to neutral. 

a * firing angle. 

X_ » commutating reactance « 2 Vfij 

O o 

D. 1 * source inductance in henry s 
s 

w direct current 
^dr " source voltage. 


(5.3) 


Rrom eqn,(3,3), 

C'^do “ •* ^dr^/ ^ ^ ^c 

*(2,34 cos 60® - 0.65)/ (| 0.1 + 0.005) 

«5.1T5P*^. 

In the derivation of eqn.C5.3) it is assumed that I^ is a 
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ripple free direct ciirreat. In the present case, the 
difference betweoi the theoretically calculated value and the 
computed value is due to the presence of ripple in the direct 
current. The ripple in the direct currsit, after allowing 
the syston to reach steady state is observed and is shown 
in fig. 3, 4, 

5'ig.3,5 shows the variation of the voltage across a 
valve over one cycle. She instants of extinction and ignition 
of valves are also shown, The voltage across valve 1 is 
indicated v/ith dotted line and the thick line variation is 
the voltage across valve 1 under ideal conditions, fhe 
ideal curve is drawn for the same overlap angle and is 
compared with the curve obtained from computer values. In 
this case also the difference observed between the two curves 
is due to the presence of ripple in the direct current. 

Pig, 3. 6 shows the variation of rectifier teminal 
Voltage and valve currents for a period of one cycle. Ihe 
valve current curve is not a dn.at topped one and the variation 
in valve current can be clearly observed which corresponds to 
the variations in dc terminal voltage, 

Ihe discrepencies in the observed waveforms and the 
theoretically calculated steady state results are due to the 
presence of ripple in the direct current. This ripple can be 
reduced by increasing the value of the smoothing reactor but 
correspondingly the settling time also increases. 
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3,4 OOFCIUSIOF 

Digital computer program is developed for the 
simulation of a 3-phase bridge converter based on the 
mathematical model given ir Chapter 2, fiis is applied, for 
the transient analysis of a converter feeding into a dc 
voltage source, She steady state waveforms of the various 
quantities such as valve voltages and currents are also 
presented and compared with theoretically calculated waveforms 
which are obtained with the assumption of ripple free direct 
current , 

The example presented in this chapter involves only 
the consideration of the normal 12 modes of the converter 
operation, The next chapter will present examples that 
involve abnormal operation of the converter. 



CHAPTER 4 


H?DO SYSTM SDroLATIOH UKDER PAIIDTED CONDITIONS 

4,1 INTRODUCTION 

The faults in DO systen can he classified as faults 
on DO line and faults in converters. The variations in 
DC link current because of faults on the DC line can be 
controlled by controlling the firing instants of valves. The 
rapid control of current is necessary as the overcurroits 
can lead to connautation failures of the inverters and damage 
to valves. Commutation failures in inverters , are generally 
selfcuring. Only under certain serious faults, the converters 
have to be blocked and bypassed. 

The faults on DC line are not considered in the present 
work but only the phenomena of commutation failure and conti- 
nuous misfire in inverter are studied, 

4, 2 SIMUMTION OP HVDG SYST®[ 

The entire HVDC power transmission systan considered 
for the present study is shown in Pig, 4,1. Por simplicity 
single 'T* network is considered to represent the tranaaission 
line. The mathematical model for the single converter used in 
the previous chapter can be duplicated for the two converters 
With the voltage source replaced by the capacitor voltage* 
Treating the two converters separately for a short interval, 
the two sets of equations cai be interlinked by the following 
equationf 




l-lVDC Pokier TRA^VSH liSl6M S/STt'^t 
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where 0^,, aad are shown in Pig. 4-1. 

She system data is taken from ref. 7 and is given below. 

AO system and transformer representation for each converter: 
Resistance per phase =0.5 ohms 

Inductance per phase = 2o mH 

AC system frequency = 50 Hz 

AO system voltage (l/N) = 42.8 kT 

Transmission line parameters: 

Equivalent series resistaice = 2.16 ohms 
Equivalent series inductance = 13 mH 
Equivalent shunt capacitance = 21.6 /iP. 

Por the present problem the value of smoothing 
reactor is taken as 1.0 H and the resistance (R^) is 

taken as 0,05 ohms. The firing angles for the valves in 
rectifier and inverter are taken as 15® and 150® respectively. 

Operation of the HVDC system is simulated with the 
following initial conditions: 

Oapacitor voltage Yo = 90 k7 
HO link current = 650 Amps. 

It is observed that the system reached steady state in 
46 cycles (0.92 sec,). This high settling time is obviously 
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due to the large value of smoothing reactor, hut to reduce the 
ripple in the direct current, the value of \ is taken as one 
henry. 

After allowing the system to reach steady state, two 
disturbances are simulated - 1. commutation failure in the 
inverter and 2, continuous misfire in the inverter. Both the 
examples are presented in the following section along with 
the results obtained, 

4.3 exahpl:es 

4,3,1 Oomrautation failure: Commutation failure is the most 
common misoperation of inverter. It is not a fault of the 
inverter but results mostly from causes that are external to 
the converter bridge. Commutation failure is the failure to 
complete commutation before the commutating emf reverses. It 
takes place because of a sudden reduction in AC voltage or 
late ignition of valves. The effect of a commutation failure 
is to reduce the bridge voltage to zero for, a short period 
during which the valve currents may increase due to discharge 
of the line capacitance. A single commutation failure is 
s elf curing particularly wh SI aided by the constant extinction 
angle controller. Cascading of commutation failures is to be 
avoided to limit the overcurrents on DC link, la double 
commutation failure, there is a time period when the inverter 
voltage reverses resulting in overcurrents on DC line. 
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In the present example a repetitive double ccmmutation 
failure is observed due to the suddm reduction in the AO 
voltage at the inverter. Under steady state conditions the 
AO voltage at the inverter end is 42.8 k7 (L/F). This is 
suddenly reduced to 15 k7 (L/F) and the observed resiilts are 
shown in I'igs, 4.2 and 4.3, with ^d without changes in valve; 
firing angle. 

Results and Discussion; Rig. 4. 2 shows the variations of 

a) Direct current at inverter terminals ^ 

b) Capacitor voltage (Y^) 

c) Inverter terminal voltage 

The purpose of the pres oit simulation is to study the 
commutation failure in inverter v;ith and vdthout the changes ; 
in valve firing angle. The results in Rig. 4.2 are observed 
when there is a reduction in the AC voltage at inverter 
station. This reduction is introduced at the end of the 
first commutation period in the 48th cycle. As a result the 
instantaneous terminal voltage is suddenly decreased from 
85.3 kY to 26.9 kY (point P.j). At this stage valve 5 and 6 
are in conducting state. At point P 2 valve 1 is fired. The 
porticai of the curve from P 2 to P^ indicates the unsuccessful 
commutation and valves 5, 6 and 1 are in conducting state* 

At the Old of this period valve 2 is fired and the terminal 
Voltage collapsed to zero as four valves (5, 6, 1 and 2) are 
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coad^ictiag. Before reaching point valve 1 ceased to 
conduct, thus a short circuit on DG link is still maintained 
hy Yalves 5 and 2, At valve 3 has to fire hut did not fire 
because of the zero voltage across it. At P^ valve 2 ceased a 
5 and 6 are in conducting state. Here the point to he noted i 
the double commutation failure (valve 5 to 1 and valve 6 to 2) 
in the inverter. Eventhough the combination 5 and 6 is a norm: 
one, as it occurred in the abnormal portion of the cycle, the 

i 

terminal voltage ^di became negative remained negative till th^ 
end of that cycle. Because of this negative voltage valve 4 
did not fire and the change in the state of the converter 
occurred only when valve 1 is fired at point Pg. Pg corres- } 
ponds to 1*2 previous cycle and the same phenomena is 

repeated. Because of no current control action and the 
discharge of capacitor whai Y^^ is zero, the direct curr^t 
increased at a fast rate and reached 2200 amps in two cycles 
aftec creating the disturbance. It is to be noted that the 
initial steady state value of the current is 630 amps, 

Pig,4.3 shows the variation of the DC quantities under 
control action at the inverter. Here the valve firing angle 
is also reduced simultaneously with the reduction in AO 
voltage at inverter station, The firing angle is reduced 
from 150® to 110® and an improved system response is 
observed. It can be noted that the occurrence of repetitive 




Ficj, 4-2 - Double^, commutation failuio due to i(,>duction iri AC voltacjc 
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double comniutatioa failure is avoided and direct current 
increased at a slower rate, 

4.3,2 Misfire in Inverter: Misfire is the failure of a valve 
to ignite. In thyristor valves which are suitably designed, 
misfire will not occur unless there is a malfunction in the 
firing control scheme. 3>uring this fault also the dc voltage 
of the rectifier bridge is reduced to zero, but for a short 
duration. This appears once in a cycle and may persist for 
many cycles depending upon whether the cause is taaporary or 
permanent. Continuous misfire of inverter valves is a serious 
fault and it impresses AG voltages on the DO line. 

Continuous misfire of the inverter valves is simulated 
by inhibiting the firing pulses to all the valves at the time 
when valves 5 and 6 are conducting. The results obtained are 
shown in Pigs. 4.4 to 4,8, 

Results and Discussion: 

Pigs. 4.4 to 4,8 show the variationsof the following 
quantities: 

a) Direct current at rectifier terminals 

b) Direct current at inverter terminals (l^l_). 

c) Capacitor voltage (V^). 

d) Inverter terminal voltage 

e) Rectifier terminal voltage - 

Prom Pigs. 4. 4 and 4,5 y it can be noted that the line 
current became oscillatory and a maximum of 8500 amps. , at 
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*■<1 ref. 


1+SI 


J'iS-4.9 Simplified Control Scheme. 

If the measured current (t < -i 

order (I ^ ^-v, ^ 

'■^d ref^> controller adranoes the ft • 

.r.w ^vances the firing of the 

valve (decreases the firing anuria') 4 -^ 

® the rectifier 

internal voltage, thus leadinxr + * 

e , nus leading to increase in I,, if the 

measured current is greater th«n +v. 

In.. « 4 .. controller 

Increases the firing ane-i. +u 1 u. 

iring angle so that the internal voltage 

decreases and hence 

I-rom the above dlsonssion, it is olear that k shonld 
negative. The follo^ng values for k and . are ohoaen 
arbitrarily , for this simplified scheme. 

k a -0,1 , 1=1,0. 

The current order ^ -•« tav.. oo 

dref baleen as 630 amps., the steady- 

state current observed from the previous examples. The 

following equation oan be easily derived from the above block « 

diagram. 


If (Ac) , tl^^^^-I^)k _ del 


“ii+1 ” “a + to 


(4.2) 

(4.5) 


where n represents the cycle number, 
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To implonent this scheme it is aeceesary that 
eqn.(4.2) is to he solved simultaneously v/ith the other 
equations defining the state of the system. In this simu- 
lation a is updated every cycle starting -with the valve 
number 2 by using the equation (4.3). 

Assuming the steady state conditions in the beginning 
the Commutation failure and continuous misfire in inverter 
are simulated and compared with the results in the previous 
section. 

Results and Discussion; 

i. Commutation failure: Pig, 4.10 shows the variation of 

a) direct current at inverter terminals (I..) 

di 

b) capacitor voltage (V ) 

c 

c) inverter terminal voltage 

The effect of control scheme seems to be negligible 
in the first two cycles and there is not much difference 
observed in the results presented in Pigs. 4.10 and 4.2. 

li. Continuous misfire t The steady state conditions 
axe taken as initial conditions and valves 4 and 5 are made to 
conduct through out the cycle, by inhibiting the firing pulses 
to other valves. results are given for 15 cycles after 

the creation of this disturbance at the end of which the 
DC link current became discontinuous at the rectifier end. 
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Pigs, 4.11 "to 4»15 show "bhe variations of 
a) Direct carreat at inverter terminals (l^^) 
h) Capacitor voltage (V^) 
c) Inverter terminal voltage 

It can he observed from Pigs. 4.11 and 4.5 that the 
inverter te3?minal current reached a maximum of 6557 amps., 
whereas it reached to 8610 amps in the same time in the 
absence of control. Goirespondingly a decrease in the line 
voltage is also observed, and can be concluded from Pigs. 4. 12 
and 4.6. The maximum voltage observed on the capacitor is 
803 kV in the presence of control and 988 kf in the absence 
of control. The current controller has little effect on 
the DC voltage at the inverter terminals. 

4. 4 COFCIUSIOF 

In this chapter HVDC system is simulated using the 
computer program developed in the previous chapter and two 
cases of misoperation of inverter valves are studied. In 
particular a study of the overvoltages and over currents on 
DC line caused by 1) commutation failure a:'d 2) continuous misfire 
in Inverter is made. In the first example a double repetitive 
commutation failirre is observed due to a sudden reduction of 
A0 voltage at inverter station.. This created overcurrents 
on DO line. These overcurrents are controlled and the astern 
performance is improved by reducing the firing ang£l© of the 
inverter valves. As second example, continuous misfire in 
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the inverter, which leahs to severe AG voltages on DC line 
is considered, This can he mitigated to some extent "using 
a current controller at the rectifier. 

An attempt has heen made to simulate the action of 
current control at the rectifier using a simplified scheme. 
The values of the gain and the time constant are chosen 
somewhat arhitrarily and it is observed '-that with this 
controller, the overvoltages and the currents are reduced 
in the case of misfire of inverter valves. The curroit 
controller has little influence on the commutation failure 
at the inverter as this event occurred within very short 
time after the disturbance (reduction in AC voltage). 



CHAPTER 5 
GOHOLUSIONS 

5.1 STMM4PX 

The dynamic simulation of the bridge con- 7 erter is 
extremely difficult because of its varying topology, caused 
by the commencement and cessation of valve conduction. ITo 
simple methods are so far available in the literature. 

In this thesis an attempt has been made to simulate 
the complete range of operation of the bridge converter on a 
unified basis. Por this purpose a detailed mathematical model I 
is developed in Chapter 2, A computer program based on -this is i 
utilized for the transient analysis of a singld converter. 

Later a HVDC power transmission system, consisting of two 
converter stations, is simulated to study the overcurrents 
and overvoltages in the system caused by commutation failure 
in the inverter and misoperation of inverter valves respec- 
tively. 

In this thesis the approach adopted for simulation is 
along the lines given in ref, 11, but the basic equations are 
derived on the basis of topological considerations. In ref. 11, 
the authors did not consider the smoothing reactor. They 
analysed only the normal operating modes and hence the 
abnormal modes corresponding to the faulted conditions cannot 
be simulated. These modes are also considered in the present 
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simulatioa, utilizing the generalized equations for the 
abnormal modes. 

5 , 2 SCOPE FOR PIIRIHER WORK 

The simulation procedure used in the present work 
can be greatly improved, if a suitable teclinique, to 
directly formulate the equations for a given mode from the 
basic equations derived in Chapter 2, is adopted. Ihe 
simulation should be in such a way that no preformulation of 
differential equations defining the converter state is needed. 

In simulating the HVPC transmission qjrstem, for 
simplicity a single T network is used to represent the trans- 
mission line. It is not a difficult task to consider several 
I sections, if greater accuracy is required or very long lines 
are considered. This will simply increase the number of 
equations to be solved* The inclusion of filter circuits and 
valve damping circuits is an extensionto the present work, but 
would not alter the basis of the program. It is to be noted 
that the representation of valves by ideal switches is not 
possible when considering the behaviour of the system at high 
frequencies. 
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